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Abstract. We discuss rank 2 sub-Riemannian structures on low¬ 
dimensional manifolds and prove that some of these structures in 
dimension 6, 7 and 8 have a maximal amount of symmetry but no 
integrals polynomial in momenta of low degrees, except for those 
coming from the Killing vector fields and the Hamiltonian, thus 
indicating non-integrability of the corresponding geodesic flows. 


Introduction 

A sub-Riemannian (SR) structure on a (connected smooth) mani¬ 
fold M consists of a completely non-holonomic (or bracket-generating) 
vector distribution A G TM and a Riemannian metric g G r(S'^A*) 
on it. For points x,y ^ M denote by 'H{x,y) the space of integral 
(horizontal) curves 7 : [ 0 , 1 ] M, 7 G A, joining x to y. 7(0) = x, 
7(1) = t/. It is nonempty by the Rashevsky-Chow theorem. 

The length functional Igi^i) = Jq ||7||gdt on the space of horizontal 
curves dehnes the sub-Riemannian distance on M by 

dg{x,y)= inf /^(y). 

'y£H(x,y) 

A curve 7 G R is called geodesic if it locally minimizes the length 
between any two (close) points with respect to dg. The description of 
most geodesics (normal ones) is given by the Euler-Lagrange variational 
principle. There is a Hamiltonian reformulation of this principle, called 
the Pontrjagin maximum principle jPMPj . It allows one to consider the 
sub-Riemannian geodesic flow as the usual Hamiltonian flow on T*M 
(abnormal extremals play no role in this respect and will be ignored 
in this paper). We will recall this together with the other relevant 
material in Section [TJ 
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As in the standard theory of Rieniannian geodesics, the metric g is 
integrable if this Hamiltonian flow is integrable on T*M in the Liou- 
ville sense, i.e. there are almost everywhere functionally independent 
integrals A = if,/2,...,/„ that Poisson-commute = 0 (see 

[Al lAKN] and also [BFj for a review of methods and problems). 

In this paper we investigate certain aspects of integrability of SR- 
structures on vector distributions of rank 2 (the smallest rank in non- 
holonomic mechanics). In general, SR-structures need not be inte¬ 
grable. For the first time, this was illustrated with a precise example 
in |MSSj by Montgomery, Shapiro and Stolin. More examples can be 
found in [Krj . We will focus on left-invariant SR-structures on Carnot 
groups, which serve as tangent cones (nilpotent approximations) for 
general SR-structures. In Riemannian geometry, the tangent cone is 
the Euclidean space and it is integrable. This integrability does not 
carry over to the sub-Riemannian case. 

We discuss integrabilit'^ of SR-structures and particularly pose the 
specihc question whether it is related to the amount of symmetry 
present in these structures. On Carnot groups of dimensions up to 
5 the geodesic flow of all left-invariant SR-structures are Liouville in¬ 
tegrable (see Section [ 2 ]), however starting from dimension 6 we show 
that the final polynomial integrals, required for Liouville integrability, 
cease to exist at least in low degrees (up to 6), even in the maximally 
symmetric situations. For precise formulations in dimension 6, 7 and 
8 , see Theorems [H [ 2 ] and [ 3 ], respectively. 

In Section [HI we reduce the corresponding systems of PDEs to sys¬ 
tems with 2 degrees of freedom in a convenient form that allows us 
to consider obstructions for integrability in a uniform setting. The re¬ 
duced systems provide a parametric 3 -components first order system of 
ODEs. Its dynamics is interesting in its own right (we speculate that 
the case corresponding to dimension 6 is similar to a forced pendulum). 

In Section [ 9 l we complement our results with the trajectory portraits 
that demonstrate irregular dynamics. Our computations show that the 
systems exhibit chaotic behavior for various values of parameters in the 
reduced formulation, providing more evidence of non-integrability. In 
dimension 8 our study agrees with the numerical observations of [Saj . 

The combination of established low-degree non-integrability, the re¬ 
duced formulation (the known integrable quadratic Hamiltonians with 

^We consider the integrals that are analytic in momenta. For a quadratic Hamil¬ 
tonian H, the existence of such an integral implies by [Dari IWh] the existence of 
an integral that is homogeneous polynomial in momenta. Moreover, in all our cases 
we need only one additional integral / commuting with H and the linear integrals, 
so this I can always be assumed homogeneous polynomial in momenta. 
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2 degrees of freedom have integrals of deg < 4), and the numerical 
evidence strongly suggests that generic SR-structures are in general 
not Liouville integrable with analytic in momenta integrals. In the 
Riemannian setting this was recently proved in [KM2] . 

The technique we use in sections [5l [6] and [7] is inherited from the work 
[KMl] , where it was exploited to prove rigorously non-existence of low- 
degree integrals for the Zipoy-Voorhees metric from general relativity 
(for related work on this topic, see |LGl IMPS! IVlj h We will explain 
the method in detail in Sections [3] and 01 In short, it allows us to 
reduce the search of integrals of a hxed degree d to a linear algebra 
problem, namely to a computation of the rank of a matrix with the size 
polynomially growing with d. The entries of this matrix are integers, 
and the computer verihcation, solely based on evaluation of the rank, 
gives a rigorous proof of the result. To the best knowledge of the 
authors, it is at present the only method that allows one to make non¬ 
existence statements for the class of integrals under consideration. 

Acknowledgment. We are grateful to the referees for useful sug¬ 
gestions improving the exposition. One of them proposed an idea of 
modular computations that we implemented in the revision. We also 
thank Vladimir Matveev for encouragement throughout our work. 

1. Hamiltonian systems and Sub-Riemannian structures 

Let us recall some basic facts from the theory of integrable systems 
and sub-Riemannian geometry. 

(A). Integrable Hamiltonian Hows. Let O) be a symplectic man¬ 

ifold. A Hamiltonian vector held is the H-dual to an exact 1-form dH : 
•) = dH. The Poisson bracket is then {F,G} = Xc). 

The Hamiltonian system is called Liouville integrable if in addi¬ 
tion to Ii = H there are integrals J 2 ,..., which are in involution, 
{Ij,Ik} = 0 , and which are functionally independent almost every¬ 
where. By the Liouville theorem a full measure set W C hP is then 
foliated by cylinders (tori in the compact case), and each cylinder has 
a neighborhood with coordinates 0 G x W, J = J{I) G M” such 
that Q, = dJ A d(j), {/ = const} ~ 'f"'"'' x M'" and the how is 0 = <h(/) 
for some function $. The noncompactness rank r can be positive for 
sub-Riemannian geodesic hows, even on compact manifolds. 

It might happen that in addition to involutive integrals Ji,...,/„ 
there are some other integrals In+i, ■ ■ ■, In+k- Then the motion is re¬ 
stricted to sub-cylinders in x - this is the resonance (the non¬ 
compact case is more subtle). Existence of additional integrals (more 
than n commuting ones) is usually referred to as super-integrability. 
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Let (M, g) be a Riemannian manifold. The geodesic flow ipt TM ^ 

tt® 

TM is the Hamiltonian flow on T*M TM with the standard sym- 
plectic structure hi if we choose the Hamiltonian H = |||p|p, p € T*M. 
The metric g is called integrable if the flow pt is Liouville-integrable. 
If M is compact, then the noncompactness rank r = 0 . 

(B) . Pontrjagin maximum principle. Consider now the non-holonomic 
case. We start with an arbitrary completely non-holonomic distribution 
A G TM and a smooth field of Riemannian metrics g^ G S^A*, x G M. 

Let {T*M, H) be the cotangent bundle equipped with the standard 
symplectic structure. The sub-Riemannian metric g yields the isomor¬ 
phism jj^ : A* —)■ A. Consider the inclusion i ■. A ^ TM. The 
following composition defines a vector bundle morphism 

T*M 4 A* 4 A 4 TM. 

Contrary to the Riemannian situation this is not an isomorphism. In¬ 
deed, we have: Ker(\['g) = Ann A and CoKer(Tg) = TM/A. 

Define the Hamiltonian function H on T*M as the composition 

T*M 4 A* 4 A ^ 4 ' M. 

This function can be locally described as follows. Let ^i,..., be some 
orthonormal basis of vector fields tangent to A. Every vector field is a 
hber-linear function on T*M. So we have H = ^ 4 - 

The Pontrjagin maximum principle |PMP] states that trajectories of 
the Hamiltonian vector field Xh in the region {iL > 0 } projected to 
M are extremals of the corresponding variational problem. They are 
called (normal) geodesics. 

Example . Consider the Heisenberg group G = Heiss with the stan¬ 
dard left-invariant metric on A = {dx^ +2^243,42) C TM.^{xi,X2,X3). 
Solving the Hamiltonian equation for 2 H = (pi -f X 2 PzY + P2 
that SR-geodesics are spirals in the direction of the x^-axis, projecting 
to arbitrary (including radius 00) circles on the plane M^(a;i,a;2). 

(C) . Vector distributions. Given a vector distribution A G TM we 
define its weak derived flag by bracketing the generating vector fields: 
Ai = A, Aj+i = [A, Aj]. The distribution is non-holonomic if A C A2 
and completely non-holonomic if A^ = TM for some k. We will assume 
that the rank of the distributions Aj is constant throughout M, then 
(dim Ai, dim A2,..., dim A^) is called the growth vector of A. 

The family of graded vector spaces {g* = Aj/Aj_i}, equipped with 
the natural bracket induced by the commutators of vector fields, form a 
sheaf of graded nilpotent Lie algebras g = 0i © • • • © over M. In this 
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paper we consider only the strongly regnlar case, when it is a bnndle 
(i.e. the strnctnre constants in the hber do not depend on a; G M). The 
typical hber is then called the Carnot algebra of A. 

For the rank 2 distribntion A C TM the prolongation is dehned as 
follows [AKl IMon] . Let M = PA = {(x, £) : x G M,i C Aj,} be the 
natnral S^-bnndle over M with projection tt : M —)■ M. Then the 
prolonged distribntion A C TM is given by C 


Example. The prolongation of the tangent bundle o/R^ is (Heiss, A). 
The prolongation of (Heiss, A) is the Engel structure discussed below. 


Even thongh the SR-behavior can be qnite different, the prolonged 
distribntion has the geometry readable off the original distribntion and, 
starting from dimension 5 , we will assnme that A is not a prolongation 
of a ra nk 2 distribntion from lower dimensions. 


2 . SR-structures on Carnot groups of dimension 3 to 5 

In this section we discnss left-invariant SR strnctnres on low-dimen¬ 
sional Carnot gronps. We will show that np to dimension 5 all of them 
are Lionville integrable. This holds for distribntions A of all ranks, bnt 
since the concern of the paper is rank(A) = 2, we restrict to this case. 

A Carnot gronp G is a graded nilpotent Lie gronp, with its Lie alge¬ 
bra 0 = fli © • • • © 0A: being bracket-generated by 0i, and distribntion 
A C TG corresponding to it. Eqnipped with a left-invariant Riemann- 
ian metric on A, snch a gronp naturally serves as a tangent cone at a 
chosen point of a general SR-strnctnre, see e.g. [BRj for details. 

In what follows we nse the following notations. A basis e* of 0 gen¬ 
erates the basis cn, G ( 0 *)* of linear fnnctions on 0*, given by 

^i{p) = {P,ei), peg*. 

We identify Cj with the left-invariant vector helds on G, and cjj with 
the left-invariant linear fnnctions on T*G. Their right-invariant analogs 
will be denoted by /* and 9 i respectively. 

The Lie-Poisson strnctnre Vlp on 0* indnces the Poisson strnctnre 
(V LP, —V Lp) on 0* © 0* and this yields the following commntation 
relation of the above fnnctions with respect to the canonical symplectic 
strnctnre on T*G: If [e*, Cj] = chcfc, then 

0 , c^j9k. 

(A). Dimension 3 : the Heisenberg SR-structure. In dimension 3 the 
only Carnot gronp@ is G = Heiss. The Carnot algebra is fieiSg = 0i©02 
with 01 = (61,62), 02 = (es) and the only relation [61,62] = 63. 


^Left-invariant SR-structures on 3D Lie groups are considered in Appendix [A] 
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The Hamiltonian H = + a;|) has two integrals: I2 = 0i and 

the Casimir = 63 = uj^. In the coordinates of Example from ^(B) 
we have: I2 = Pi, I3 = Ps- These Ji = H,l2,l3 are involutive and 
functionally independent, and so yield Liouville integrability. 

There is also a fourth (noncommuting) integral 14 = ^2 = P2 + Xip^ 
conhning the motion to the cylinders x C T*G = G x g*, and 
making the system super-integrable. 

Actually, for all systems considered in this paper whenever we es¬ 
tablish Liouville integrability, the super-integrability (but not maximal 
super-integrability) will follow. Indeed, we will always indicate a right- 
invariant linear form (commuting with the left-invariant Hamiltonian) 
that is functionally independent of the other integrals. 

(B). Dimension 4: the Engel SR-structure. In dimension 4 we also 
have only one SR-structure, related to the well-known Engel structure. 

The graded nilpotent Lie algebra is 0 = 0i © 02 © 03 = (61,62) © 
(63) © (64) with the nontrivial commutators: 


[61,62] — 63, [61,63] — 64. 

The Hamiltonian is H = + wl), and h = O2, I3 = 63, R = 9 ^ 

together with R = H form a complete set of integrals. Adding R = 61 
makes the Hamiltonian system super-integrable (notice though that R 
does not commute with RyR)- 
In coordinates on G we hav^: 

2H = (pi + X2P3 + X3P4f + pI, 
and the integrals are: 

R = 62= P2 + X1P3 + Ixjpi, R = e^= P4; 

R = 93 = p 3 + X1P4, {R = - 9 i=pi). 

Alternatively, to get an involutive set of integrals, we can use the inte¬ 
grals J2 = R, J3 = R and the Casimir function J4 = — 2u:2u:4\ 

J2 = Pi, J3 = Pi, Ja= p\ — 2 p 2 P 4 = — 2/2/4. 

The obtained integrals establish Liouville integrability of H. 


^This and similar formulae are obtained via realization of the basis as left- 
invariant vector fields on G. For the Engel structure: ei = —{dx^ +X 2 dx,i +X 3 dxi), 
62 = 9x2, 63 = 9x2, 64 = 9xi- The right invariant vector fields fi are such fields on 
G that commute with Cj and have the same values at the unity of G. 
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(C). Dimension 5 : the Cartan SR-structure. In dimension 5 there are 
two SR-structures: one on the prolongation of the Engel structure (a 
partial case of the Goursat structure, easily seen to be integrable, so we 
skip it) and the other related to Cartan’s famous (2, 3 , 5 ) distribution. 
The Carnot algebra is the positive part of the grading, corresponding 
to the first parabolic subalgebra of the exceptional Lie algebra Lie(G2): 
0 = 01 ©02 ® 03 = (ei, 62) © (es) © (64, 65). The nontrivial commutators 
of 0 are: 

[61,62] = 63, [61,63] = 64, [62,63] = 65. (1) 

The Hamiltonian is if = + ^2)- terms of right-invariant vector 

helds and the corresponding linear functions, we have the following 
involutive set of integrals: the Casimir function 

I2 — ^1^5 ~ ^2^4 + 1^3 = ~ UJ 2 ^ 4 : + ^^31 

together with the linear integrals I3 = 6*3, J4 = 6*4, R = 6*5. Again 
adding Jg = 61 makes the Hamiltonian system super-integrable (the 
next obvious candidate Jg = 6*2 is already functionally dependent with 
the previous integrals; they do not commute with I2). 

In coordinates on G we have: 

2if = (pi - lx 2 P 3 - XiX2PAf + {P2 + \xxP3 + XiX 2 P 5 )^, 
and with the notation J± = Xip^ ± X2P3 the integrals are: 

h = P1P5 - P2P4 + |f»3 + + IP3J+, h = Pa, 

h = P 3 , h = P 5 - 

The additional integral is either Jg = Pi + ^X2P3 + {x3 — ^XiX2)pA + ^X2P5 
or r^=p2- ^xips - \xIpa + (0:3 + \xiX2)p3- 

3. Discussion: on detecting integrability 

As discussed in the previous section, all left-invariant SR-structures 
on Carnot groups of dimensions < 6 are Liouville integrable, and 6 is 
the smallest nontrivial dimension from this viewpoint. Reference |MSS] 
provides an example of a non-integrable left-invariant SR-structures on 
a Carnot group in 6D. Namely, it is supported on a rank 3 distribution 
with growth vector ( 3 , 5 , 6 ) on the flag manifold SL( 4 )/R, where B 
is the Borel subgroup. In the respective | 3 |-grading st( 4 ) = 
with b = XlLofl* distribution A corresponds to 0_i. Integrability 
in loc.cit. is understood in algebraic sense, and the obstruction to it is 
based on a result of Ziglin (using the monodromy group). 

We study SR-structures on rank 2 distributions of other 6D Carnot 
groups G, and especially ask whether integrability is linked to the max¬ 
imal symmetry of a structure. Of course, every symmetry field yields. 
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by the Noether theorem, an integral linear in momenta (to be called, in 
what follows, a Noether integral), bnt the qnestion is if these (plentifnl 
in onr cases) linear integrals snffice for integrability. 

Here, we consider two kinds of inhnitesimal symmetries. On the 
one hand, we consider vector fields whose flow preserves the nnderly- 
ing rank 2 distribntion of the SR-strnctnre. On the other hand, we 
consider Killing symmetries whose flow, in addition, preserves the SR- 
metric. In both cases we observe that the existence of an additional 
symmetry appears not to be linked to the existence of a complete set 
of involntive integrals. In each of onr examples, considering a sys¬ 
tem with D = dim(G) degrees of freedom, we will be able to identify 
D — 1 involntive integrals (1 Hamiltonian and D — 2 Noether integrals), 
henceforth concentrating on the search for the final integral. 

In dimension D = 6, we will list all left-invariant SR-strnctnres on 
Carnot gronps and discnss their integrability. It will be shown that a 
hnal integral of low degree does not exists even if A is the maximally 
symmetric distribntion; on the other hand for the maximally symmetric 
SR-strnctnre snch an integral does exist. In dimension D = 7 or 8 we 
will focns on maximally-symmetric SR-strnctnres. These have D + 
1 linear independent Killing symmetries (the corresponding Noether 
integrals are not in involntion). This follows from the Tanaka theory 
that we review in Appendix |Bl We can also straightforwardly hnd the 
Noether integrals by realizing the gronp law via the Baker-Campbell- 
Hansdorff formula and solving the relevant system of PDEs. 

Non-integrability, on the other hand, is more difficult to prove. There 
are few methods to detect it, depending on the integrability setup. An¬ 
alytic non-integrability on a compact manifold follows from positivity 
of the topological entropy, see [T] (this is guaranteed if a horseshoe 
subsystem exists, which can be obtained via scattering data, see [B]). 
but this cannot be applied in our case. 

Obstructions for algebraic integrability can be found via differential 
Galois theory |MR] . and the Painleve test. The latter assumes the inte¬ 
grals to be rational in all coordinates. In contrast, we are interested in 
first integrals that are smooth in the base coordinates and polynomial 
in momenta (such integrals, also known as Killing tensors, have appli¬ 
cations in mathematical physics and general relativity), and the above 
tests are not applicable. The method that can detect existence of this 
type of integrals was proposed in the work |KMlj . Before describing it 
in details in the next section, let us explain the simple idea behind it. 

The condition governing existence of an integral of degree d is an 
overdetermined system of linear differential equations on 
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unknown functions of D variables. It is of finite type, meaning the 
system is reducible to ODEs. 

Checking the explicit compatibility conditions can be cumbersome. 
Instead, we compute all differential consequences, i.e. take partial deriva¬ 
tives of the PDEs in the system, and it is enough to take those of total 
order < d + 1 , cf. |Wol IKMlj . At this stage, the problem of solving 
differential equations is reduced to linear algebra: kernel of the matrix 
of the system, evaluated at a hxed point, corresponds to degree d inte¬ 
grals. Some of those are products of apriori known linear and quadratic 
integrals, which gives the lower bound on the nullity of the matrix. 

The decision on existence of the hnal integral is thus reduced to the 
computation of the ra nk of this matrix. If the nullity is the minimal 
possible (decided by the rank), no additional integral exists. 

Similar to applications of the Galois theory or Painleve test, our 
method can be implemented on a computer. Only integer numbers are 
involved, so the symbolic computations are exact, containing no round¬ 
ings or approximations (with large numbers involved, these computa¬ 
tions are more reliable than a calculation by a human). Signihcance of 
such proofs in mathematics has steadily increased in recent years. 

In the next section, we give the technical details of the method, and 
then use it to prove Theorems [H [ 2 ] and [ 3 ] in a mathematically rigor¬ 
ous, computer-assisted manner. The reader interested in independent 
verihcation of our computations is invited to download the Maple-code 
from the supplement to arXiv: 1507 . 03082 v 2 . 

4. Method to check existence of the final integral 

Similarly to |KMlj . we show that the systems, considered further in 
this paper, lack one hnal integral F required for Liouville integrability. 
In all our cases, with D = dim G degrees of freedom, we have D — 2 
commuting Noether integrals that we normalize to be ps, ■ ■ ■ ,Pd- 

This reduces the system to a (D — 2 )-parametric Hamiltonian system 
with 2 degrees of freedoms (simply let pi = const*). It means that the 
unknowns (=coefhcients of F) depend on 2 variables only. 

We however do not perform the reduction until Section [HI and keep 
the momenta p = {pi,... ,pd) in the general form. The Hamiltonian 
thus takes the form H = H{xi,X2,Pi, ■ ■ ■ ,Pd)- 

A hrst integral that is smooth by the base variables x = {xi,X2) 
and polynomial in momenta of degree d, and that commutes with the 
Noether integrals P3,... ,Pd, has the form 



(2) 
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(jf = nliiP? for a multi-index r = (ri,..., td) and |r| = '^0- 

The requirement that functions of the form F are integrals is given 
by the Poisson bracket relation {if, F} = 0 , which is a homogeneous 
polynomial in momenta of degree d-l- 1 . Thus we have first order 

PDE (coefficients of {H,F}) on unknowns (coefficients Qr of 

F). Let us call this linear PDE system Sd- 


(A). The bounds on the number of integrals. Instead of treating the 
differential system Sd, we consider the associated system of linear equa¬ 
tions, given by fixing a point o G G. Namely, denote by the f-th 
prolongation of Sd- This is the system obtained by differentiating the 
PDEs from Sd by xi,X2 up of total order < k. The total number of 
the equations is consequently md^k = • 

The unknowns are now the derivatives a°.^ = d„ar{o) (jets) with 
multi-indices r = {ti, ... ,td) and a = {ai,a2) (of different size) of 
lengths |r| = d and |ct| < k. The collection of the unknowns is denoted 
V = and represented as a column-vector. Their number, denoted 
#1/ (the height of this vector), is equal to Ud^k = 
ffc) 

The system evaluated at o G G has the form M ■ V = b with 

ik) 

some md^k x matrix M = . Let be the number of linearly 

independent first integrals of degree d. We obviously have the upper 
bound, in which the right hand side stabilizes for k = d + 1 (cf. |Wo] I: 

Ad<^r:=#foi'^-rankMf. ( 3 ) 

Due to this stabilization let us denote in what follows: 6 d = 

Vd = and Md = 

On the other hand, our system possesses integrals G = H (quadratic: 
di = 1 ) and I2 = Ps, ■ ■ ■, Id-i = Pd (linear: dj = 1 , 1 < i < D). The 
derived integrals ]{[ J™* of degree Y = d will be called trivial. 
Thus we deduce the lower bound 


A. > AO : = 



d — 2 i + D 
D -3 



( 4 ) 


If we show that the bounds in ([ 3 ]) and ([ 1 ]) coincide, A^ = 6 d, then 
we conclude that A^ = A^, and so all integrals of degree d are trivial 
(reduced to already established Noether integrals). This gives non¬ 
existence of the final integral in degree d. 


(B). The procedure. There are two important differences to |KMlj 
that facilitate our computation. First, our model is homogeneous, so 
the choice of point is not essential (in the general case we have to choose 
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a generic point, which gives the stable values of Sd,k)- We always choose 
o = ( 0 , 0 ) in the plane Second, our Hamiltonian H (scaled 

by some integer factor that is at most 288 ) is a polynomial with integer 
coefficients, so we do not need to handle rational expressions. 

Most complications are related to the calculation of rank(M). The 
large size of the matrix makes Gaussian elimination costly (the dimen¬ 
sion of the configuration space in our computations D = 6 , 7,8 exceeds 
dim = 4 of [KMlj b But the matrix contains many zeros, and further 
simplifications are possible as follows: 

1 . In all equations all coefficients are kept integers at all stages 
by multiplying by their common denominator. Proportional equations 
(rows of M) are removed. 

2 . We have the freedom to add constant multiples of trivial integrals 
to F defined in ([2]). At the point o the corresponding unknowns can 
therefore be cleared from the equations. Let V^pA C V embrace all (su¬ 
perfluent) unknowns that do not appear in the equations. We remove 
the corresponding zero columns from M. We have: #Kpfl ^ ^d- 

3 . Perform a partial solution of the system, iteratively solving the 
monomial and bimonomial equations until no more such equations re¬ 
main. Let Kaon and Vbimon be the corresponding unknowns. 

Denote the matrix of the reduced system (obtained from Md by the 
above simplifications) by M^ed, and let Hred = Vd\ (HspsU VmonUl4imon)- 
The reduced system then reads M^ed ■ fred = 0 , and we get the formula 

Sd = Wd - rk(Md) = {Wred + #VAon + #Wimon + #Kpfi) 

(rk(Afi.ed) T ifFraon T TT^Lbimon) = 4 f^red T TT^L^pfl rk(Afi.ed)- 
Notice that #V)ed = rk(Mi.ed) and #Hspfl = Aj) imply 5 d = A^J. 

(C). The modular approach. Our algorithm confirms non-existence of 
the final integral of degree d when Sd = A^. The right hand side is 
given by (jl]), while the left hand side depends on rankM^ as in ([ 3 ]). 

To reduce the rank computation we work modulo p for a prime p. De¬ 
note by L[p] the matrix obtained from a matrix L by passing to mod p 
entries. Since for a square matrix L we have det(L[p]) = (det L) mod p, 
we conclude that 

rank Md [p] < rank Mrf. ( 5 ) 

For some specific values of p the rank can actually decrease upon com¬ 
puting modular, but for sufficiently large primes p we have equality in 
(] 5 j). Thus, if for some prime p 

^d\p\-= Wd - rankMd[p] = A(), 

then we conclude non-existence of the final integral of degree d. 
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The main complication is however to hnd such p. Our experiments 
show that the decisive p grows fast with D. So even though these mod¬ 
ular computations are cheaper for every particular p, going successively 
by increasing primes actually increases the computation time. Yet, a 
choice of a random increasing sequence of p turns out to be useful. 


5. Left-invariant SR-structures in dimension 6 


In this section we show a certain type of non-integrability for a rank 2 
left-invariant distribution on a 6D Carnot group G. Every such 2 - 
distribution A is encoded as the space g_i in the corresponding graded 
nilpotent Lie algebra g. 

In 6D the growth vector is ( 2 , 3 , 5 , 6) (recall we assumed that A is 
not a prolongation of another rank 2 distribution), and every such Lie 
algebra g is a central ID extension of the Cartan algebra from Section 
|!!](C), the distribution also being an integrable extension |AK] . 

Thus g = gi © g2 © gs ® 04 = (ei, 62) © (cs) © (64,65) © (ee) has Erst 
commutators as in (HD, which should be accompanied by the brackets 
01 © 03 —t g4. This leads to precisely three algebras, called elliptic, 
parabolic and hyperbolicQ in |AKj . We will study them in turn. 


(A). Integrability of the maximally symmetric elliptic SR-structure. 
The elliptic ( 2 , 3 , 5 , 6 )-distribution has the following structure equations: 

[61,62] = 63, [64,63] = 64, [62,63] = 65, [64,64] = 66, [62,65] = 60. 

Its symmetry algebra has dimension 8 EH. and it is not maximally 
symmetric as a 2-distribution, but it supports the maximally symmet¬ 
ric SR-structure. Namely, dehning the SR structure by the orthonor¬ 
mal frame 64,62, we conclude that its symmetry dimension is 7 (see 
Appendix IHD. The corresponding Hamiltonian is 

2H = {pi-lx2P3-XiX2P4-lxlx2PeY+{p2 + lxiPs+XiX2P5+lxixlpey. 

There are two Casimir functions R = uq and C = —u^Uq. For 

the maximally symmetric Hamiltonian H = + 0J2), these together 

with J4 = H, I2 = 0J1UJ3 — UJ2UJ4 + itUg and the right-invariant linear 
functions I3 = 63, I4 = 64, I5 = 6*5 and Iq form 6 involutive integrals 
(C = |(/| + /|) — hh), so this system is Liouville integrable. Notice 
that I2 = 6^46*5 — 9204 + is also an integral, and I2 — I2 = h ■ K, 
where K is the last Killing vector held (neither I2 nor K commute with 
Ii,... ,Iq, but they make the system super-integrable). 

^The (2,3,5,6)-distributions are given by a conformal quadric on gi due to con¬ 
formal identification adgj : 0 i ~ 33 , whence elliptic, parabolic and hyperbolic. 
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In coordinates: J3 = ps, I4 = p^, J5 = ps, Iq = pe and 

h = {Pl - \X2P3 - XiX2Pi - lxjx 2 P 6 ){p 5 + X 2 Pe) 

- {P2 + ^XiPs + X1X2P5 + \xixlp(i){p 4 : + XiPe) 

+ |(P3 + XiPa + X2P^ + ^{xl + xDpqY- 

(B). Non-integrability of the parabolic SR-structure. The parabolic 
( 2 , 3 , 5 , 6 )-distribution is given by the structure equations: 

[61,62] = 63, [61,63] = 64, [62,63] = 65, [61,64] = 66- 

With its ll-dimensional symmetry algebra it is the maximally sym¬ 
metric non-holonomic rank 2 distribution in 6D, see [DZlIlK]. 

Up to equivalence there is only one left-invariant SR-structure (this 
follows from the fact that the Tanaka prolongation g of the Carnot 
algebra g has go C gl(gi) equal to the Borel subalgebra, and the cor¬ 
responding group transforms the invariant SR-structures), and it is 
given by the orthonormal frame 61, 62 (the symmetry dimension of this 
SR-structure is 6, and so it is not maximally symmetric). The corre¬ 
sponding Hamiltonian H = + a;|) has the coordinate form 

2 H = {pi - \x2P^ - X1X2PA - ^xlx2Pey + {p2 + Ixips + X1X2P5Y. (6) 

There are two Casimir functions = 6*5, cug = 6*6, and two additional 
Noether integrals 6*3 = ps, 6*4 = p4, that form an involutive family 
I2 = P3,h = P4:,h = P5,h = Pe- However no other Casimirs or 
commuting linear integrals exist. 

In search of more complicated integrals we perform the computations 
for the hnal (6th) integral of degree d and arrive at the following result. 

Theorem 1. The final integral of degree d < 6 for the Hamiltonian ^ 
of the left-invariant SR-structure on the parabolic ( 2 , 3 , 5 , 6 )-distrihution 
does not exist. 

Proof. First let us notice that it is enough to prove non-existence of a 
nontrivial integral Iq of degree 6. Indeed, if a nontrivial integral I of 
degree d < 6 exists, then I ■ is a non-trivial integral of degree 6. 

Therefore we shall apply the procedure described in Section 0 ] to 
our system with d = 6 onljo. For sextic integrals, seven prolongations 
need to be performed in order to achieve equality for the upper bound 
6q = . Our computation gives: 


# all eqns 

#^6 

# eqns 

Wred 

rk(Mred) 

^6 

28512 

20790 

11816 

9155 

9155 

130 


^In fact, we run the test for 1 < d < 5 as well, confirming the same result. 
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The last number coincides with the number of trivial integrals 
Ag = 130 , and hence by the discussion in § 5 ] there is no integral of 
degree 6, which is independent of and commuting with I2, ... ,15. □ 

(C). Hyperbolic and other elliptic SR-structures in 6 D. The hyper¬ 
bolic rank 2 distribution with growth vector ( 2 , 3 , 5 , 6) has the following 
structure equations: 

[61,62] = 63, [61,63] = 64, [62,63] = 65, [61,65] = 66, [62,64] = 66- 

The Hamiltonian corresponding to orthonormal frame 61,62 is 

2 H = {pi-^X 2 P 3 -XiX 2 P 4 -jXixlpey+{p 2 + lxiP 3 +XiX 2 P 5 +lxlx 2 Pey. 

There are two Casimir functions R = ojq and C = uJiUJ3 — uj3UjQ. For the 
Hamiltonian H = -|- 1X2) these two together with R = H and the 

integrals I3 = O3, I4 = 64, R = 6*5 (6*6 = Ci;6) form 6 involutive integrals, 
but they are functionally dependent {C = /4J5 — I3R). 

The most general left-invariant SR-structure on both the elliptic 
and the hyperbolic ( 2 , 3 , 5 , 6 )-distributions can be brought into the form 
2 H = cjf + {aoji + 6^2)^, 6 7^ 0 , with the same 4 Noether integrals. 
However no other Casimirs or commuting linear integrals exist. 

In all these cases (except the elliptic case with a = 0 , 6 = 1 ) the 
system is neither SR-maximally symmetric (the symmetry algebra has 
dim = 6), nor maximally symmetric as a distribution (the symmetry 
algebra has dim = 8). 

In all these cases the search for the hnal integral reduces to the same 
problem. We can apply the machinery used in Theorem [H and the 
computations show the same non-existence result (in all cases except 
the elliptic a = 0 , 6 = 1 ). This non-existence of low degree integrals 
suggests that these Hamiltonians are not integrable. 

6. Maximally symmetric SR-structures in dimension 7 

A non-integrability effect established in the previous section happens 
also in higher dimensions. We noted that the parabolic distribution 
A in 6D is maximally symmetric, but for the left-invariant parabolic 
SR-structure (A, g) in 6D the symmetry algebra of (A, g) is minimal 
possible: the algebra of left-translations g. 

In general, the symmetry algebra of a left-invariant SR-structure 
{A,g) on a Carnot group G is a graded Lie algebra g and it contains 
the Lie algebra of G, namely g = gi © ... gjy C g. The additional part 
is contained at most in the zero grading: g/g = go |Moj . Clearly this 
piece is at most 1-dimensional go C so(gi,(7). 


®We provide a simple proof of this fact in Appendix [BJ 
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Thus dim Sym(A, ^f) < dimg + 1 . The equality is attained if the 
rotation endomorphism 0 G so(gi,5f) extends (uniquely) to a grading 
preserving derivation of g. Let us investigate if such a maximally sym¬ 
metric left-invariant SR-structure on a Carnot group is integrable. 

In 6D the only maximally symmetric SR-structure is the (unique 
up to scale) SR-structure on the elliptic ( 2 , 3 , 5 , 6 )-distribution (with 
dim Sym = 7 ) and it is integrable. Consider the case dim G = 7 . 

Here the only maximally symmetric SR-structure on a rank 2 dis¬ 
tribution A on a 7 D Carnot group G (that is not a prolongation from 
lower dimensions) with dim Sym(A, ^f) = 8 has growth vector ( 2 , 3 , 5 , 7 ) 
and the following structure equation^ of the graded nilpotent Lie al¬ 
gebra g = Lie(G): 

[61,62] = 63, [61,63] = 64, [62,63] = 65, 

( 7 ) 

[61,64] = —[62,65] = 65, [61,65] = [62,64] = 67. 

Here H = S'lid go = (62®a;i — 6i(8)a;2)- There are 3 Casimir 

functions cjg, C07 and U3{uI-\-Uj) — ^{uI—uI)uq—U4U5Ut. The involutive 
family of integrals 63,... ,6i generates these Casimirs and together with 
the Hamiltonian they lack 1 more integral for Liouville integrability. In 
local coordinates, we have 

2 H = [pi - \x2p3 - X1X2PA - \x\x2p% - \x1xlp7f 

+ (p2 + \x1p3 + X1X2P5 - \xixlpe \xIx2 P7Y , (8) 

and the integrals are I2 = P3, ..., Iq = P7. Looking for the hnal integral 
/7, we again invoke the method of Section 0 ] to obtain: 

Theorem 2 . The final integral of degree d < 6 for the Hamiltonian ^ 
of the left-invariant SR-structure on the (2, 3 , 5 , 7 )-distribution given by 
o does not exist. 

Proof. We perform the same computations as in the proof of Theo¬ 
rem [H In this case, our computer capacities allow us to study polyno¬ 
mial integrals up to degree d = 5 . We need six prolongations to arrive 
at a dehnite conclusion, which is presented in the table: 


# all eqns 

#^5 

# eqns Mred 

#Rred 

rk(Mred) 

^5 

25872 

16632 

9397 

6993 

6993 

166 


Since the number ^5 = 6 ^'^ coincides with the number of trivial integrals 
Ag = 166 , we conclude absence of the hnal integral of degree d < 5 . 

'^These are obtained from the (2,3,5,6) parabolic distribution by the central ex¬ 
tension technique of |AK) . 
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To handle the case of degree d = 6, we use the modular approach, 
described in Section [ ^(C) The computation concludes faster, but to 
reach a dehnite answer we need a suitably large prime. In our case 
p = 101 suffices and we obtain the following result: 


ff all eqns 

#^6 

# eqns ilA-ed 

#Ked 

rk(Mred) 

56 101 

61776 

41580 

19137 

15848 

15848 

296 


This computation implies Sd\p] = A°, and we conclude non-existence 
of the hnal integral of degree d < 6. □ 

Remark. The indicated p for d = 6 is not claimed to be the minimal 
possible. But search for the minimal p requires more computer time. 
For instance, with d = 5 the computation for d = 5 gives 6d[p] > 
for the primes p = 2,3,..., 29, and we obtain equality (implying non¬ 
existence of degree 5 integral) for the next primes p = 31, 37 and 41. 


7. On integrability oe SR-structures in dimension 8 


There are two SR-structures p on a rank 2 distribution A on a 8D 
Carnot group G (that is not a prolongation from lower dimensions) 
with dim Sym(A, g) = 9: one with the growth vector (2,3,5,6,8) and the 
other with the growth vector (2,3,5,8). The distributions are obtained 
by central extension from 7D as in |AKj . and we take the (unique up 
to scale) so(2)-symmetric metric g (in the cases, when go so (2)). 

The second SR-structure (A, g) has a more symmetric underlying 
distribution (with the symmetry dimension 12 vs. 10), but it is the 
first one that is integrable. 


(A). The ( 2 , 3 , 5 ,6,8) SR-structure. The structure equations of the al¬ 
gebra g = Lie(G) = gi©- • -©gs = (ei, e2)©(e3)©(e4, e5)©(e6)©(e7, es) 
are the following: 

[61,62] = 63, [61,2,63] = 64,5, [61,64] = [62,65] = 66, 

[61,2, 65] = 67,8, [63, 64,5] = i 67,8, 

where we use complex notations ea,b = ^a+i Cb- In this form it is obvious 
that the action of SO( 2 ) on g, composed of the standard action on gi, 
g3, gs and the trivial action on g2, g4, is an automorphism. 

The left-invariant Hamiltonian H = + U2) has 5 commuting 

right-invariant Killing helds (integrals) I2 = 0 ^, I3 = 9 ^, J4 = 9 q, 
G = 6t, Iq = 6s. In addition, there are 2 Casimir functions 

cu| + , n n ^4 F ^5 

Ij = UJiUs — UJ 2 UJ 7 + UJsUJq -^- = 6i6s — 6267 + 6 s 9 q - - —, 

C = UJ4UJ7 + a; 5 a ;8 — = 6467 + — \ 9 \, 
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of which the second is dependent on / 2 ,..., Jy. Yet we have one more 
quadratic integral 

Jg = CJiCUs — CJ2^4 + ^^3, 

and it is straightforward to check that the involutive integrals Ji = 
are functionally independent almost everywhere on T*G. 
Consequently, the considered SR-structure is Liouville integrable. No¬ 
tice that Jg = 0105 —02 6 ^ 4 - 1 -is different from Jg and is also an integral 
of H, which again manifests super-integrability. 

In coordinates, denoting + x^, we have 

UJl=Pl- \X2P3 - X1X2P4 - \xlx 2 P (3 - + xl)x 3 P^ \X1X2X3P3, 

ui2=P2 + ^XiPs + X1X2P3 + Ixixjpe IX1X2X3P7 -l{xl + a‘^)x 3 Ps, 
1^3 = P3 + XiPi + X 2 P 5 + ^P6 + (fo^l + X2X3)P7 + {^X2 - XiX3)ps, 
U 4 =P 4 + XiPq -h ^xIpj {lxiX2 - X3)ps, 

(X 5 =P 5 + X2P& + (|a;ia;2 + X3)p7 -f \xlp3, 

We = P 6 + XiP 7 + X2P8, UJj = P 7 , Wg = P 8 , 
and 6 i = pi for 4 < i < 8 ; the formulae for Jj follow. 

(B). The (2,3,5,8) SR-structure. The free truncated graded nilpotent 
Lie algebra fl = gi © • • • © fl 4 = (ei, 62 ) © (eg) © ( 64 , 65 ) © (ee, ey, eg) 
with the structure equations 

[ei,e2] = eg, [ei,e3] = e4, [e2,e3] = eg, 

[ei,e4] = ee, [ei,e5] = [e2,e4] = ey, [e2,e5] = eg 

was also studied in |Sa] . The left-invariant Hamiltonian H = + 

w|) has 6 commuting right-invariant Killing helds, leading to Noether 
integrals I 2 = O 3,13 = 64,14 = 05 , J 5 = 06, = 07 , /y = 0g. In addition, 

there is 1 cubic Casimir function, but it depends on the linear integrals. 

Thus we again lack one hnal integral for integrability. To set up its 
computation we write the Hamiltonian in local coordinates: 

2H = (pi - ^X2P3 -l{xl+ xl)p5 - \xixIp 7 - ^xlps)^ 

+ {P2 + Ixips + \{xl+ xl)p4 + \x\p(3 + \xlx 2 P 7 Y ■ ( 9 ) 

Theorem 3. The final integral of degree d < 5 for the Hamiltonian 
m of the left-invariant SR-structure on the (2,3,5,8)-distribution does 
not exist. 

Proof. We use again the procedure from SectionlUto show non-existence 
of a non-trivial integral of degree 5. After six prolongations of the PDE 
system, we arrive at the following table: 
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^ all eqns 

#1^5 

# eqns M^ed 

#V7ed 

rk(Mred) 

^5 

48048 

28512 

4439 

3514 

3514 

314 


The upper bound ^5 = 314 = realizes the number of trivial integrals, 
and proves that no hnal ( 8 th) integral of degree d = 5 exists. □ 


8 . Reduction to the system with 2 degrees oe freedom 


In this section, we give a uniform description of several of the pre¬ 
viously discussed systems in terms of first order ODE systems in 3D. 
In particular, we reformulate in this way the three systems exhibiting 
non-integrable behavior, namely the ( 2 ,3,5, 6 ) parabolic, ( 2 ,3,5,7) ellip¬ 
tic and (2,3,5, 8 ) free truncated SR-structures given by the Hamiltoni¬ 
ans (j 6 ]), ([ 8 ]) and (jO]). In addition, the same reduction can be performed 
for the general (2,3,5, 6 ) elliptic and hyperbolic SR-structures. 

First, note that in all these cases the Hamiltonian is a sum of two 
squares and so can be expressed as 

2H = cos^z / sin^z, (10) 


and Pi = Ci for i = 3,... ,D are the Noether integrals. Symplectic 
reduction via these integrals (hxing them and forgetting about Xi, 3 < 
i < D, of which nothing depends) is a classical procedure, see [WhllA]. 
Thus, in view of (ITOD . Hamilton’s equations can be rewritten in terms 
of X, y, z and p (as well as C3,..., cg). 

For instance, in the case of the parabolic (2,3,5,6)-problem, we ex¬ 
press the coordinates pi,p 2 in terms of the coordinate z as follows: 


Pi = p cos Z + IX2C3 + X1X2CA + \x\x2CQ, 

P 2 = psinz - \x 1 C 3 - X 1 X 2 C 5 . 


Next, we can conhne to an energy shell, that is £x R = = const, 

which reduces the dynamics to the manifold S'iT*M^ = R^(a;, y) x 5 ^( 2 ;), 
where we let x = xi, y = X 2 . Without loss of generality we can assume 
p = 1. After an appropriate change of coordinates, the Hamiltonian 
equation 7 ) = {p, H} on the energy shell writes as the 3x3 system: 


X = cosz, y = sin z, z = Q{x,y), 


( 11 ) 


where Q = Q{x, y) is a quadratic polynomial. This polynomial can be 
brought to the following normal form (a 7 ^ 0 & 6 7 ^ 0 ) 

Q = Qi{x^y) = ax^ -\-hy for D = 6 parabolic, 

Q = Q 2 {x,y) = ax^ + by"^ + c for D = 7, 8 


( 12 ) 

(13) 
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(the latter formula contains also the 6 D elliptic and hyperbolic cases). 
Take, for instance, the 6D parabolic case, formula fll2l) . In this exam¬ 
ple, we have a = Cq/2 and b = C 5 , and we assume a, 6 7 ^ 0 . 

Notice that the condition a, 6 7 ^ 0 is important, as otherwise the 
system hbers over a 2D flow, which can never be chaotic. 

A similar effect happens for a = b and Q = Q 2 {x, y), where a change 
of variables x = r cos'ijj, y = rsin'ip reduces the system to a 2D flow 
with coordinates r and s = z — 'ip- This latter case corresponds to the 
6 D elliptic maximally symmetric SR-structure. The corresponding 3D 
system possesses the following integral 

F = -|- I — rsin(s), 

which corresponds exactly to the integral I 2 identihed in Section |5l cf. 
also na. However, for the general a, b, it will be shown in the next 
section that the system exhibits a chaotic behavior. 

Remark. One can check that in the complement to a hypersurface the 
following 1-form on Mf{x,y) x is contact: 

a = ^{ax^dy — b y^dx) + ^{xdy — y dx) + cos zdx + sin 2 ; dy. 

In this domain its Reeh vector field Ra, given by the two conditions 
a{Ra) = 1 and da{Ra, •) = 0, preserves the volume form a A da and 
so is divergence-free with respect to it (the Reeb field Ra plays a distin¬ 
guished role in contact geometry). Our vector field, given by HU) for 
Q = Q 2 {x,y), is proportional to Ra, and so has the same trajectories. 
For Q = Qi{x,y) the situation is similar, if a is properly modified. 

We conclude this section with a note on the resemblance of system 
(ITT]) to a driven pendulum in the case Q = Qi{x,y). Let us eliminate 
x,y from (ITT]) . Differentiating i and replacing x and y via ODE (ITT]) , 
we get the following 3rd order ODE on 2 : = z{t), where A = ^ o 

A ( 2 ;" — b sin z) = 2a cos 2 ;, 

which can be written in non-local form as: 

z" — b sin 2 ; = A“^ (2a cos z) = 2a cos 2 ; cos 2 ;. (14) 

In this form it resembles the driven pendulum z" — bsinz = a cos kt 
without dissipation. For a = 0 system (ITT)) is the simple pendulum 
when b < 0, while for b > 0 the second term on the left hand side de¬ 
scribes a repulsive poweiH. However, contrary to the driven pendulum, 
where the right hand side is an external force, system fll4p seems to be 
self-driven. The evolution of this system is shown in Fig. [T]for three 
different parameter combinations. The orbital dynamics in Fig. [T] is 

®For instance, when z -C 1 the solutions are hyperbolic. 
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quite complex and resembles the dynamics of the damped driven pen- 
dnlum (see, e.g., Fig. 9, 10 in m and references therein), indicating 
non-integrability. This resemblance appeals for a more systematic nn- 
merical analysis of system (ITT]) , which is provided in the next section. 


*N 



Figure 1. The orbital evolution of variables z, z' for the 
D = 6 parabolic case with the parameters a = 10 and 
b = —0.1 (left panel), b = —1 (middle panel), 6 = 1 (right 
panel). The initial conditions for {x, y, z) are (0, —5, 0) in 
the left and middle panels and (0, 0, 0) in the right panel. 

The red curves show the evolution in the time interval 
0 < t < 1, while the black dots continne it to time 1000. 

9. Numerical evidence of non-integrability 

In this section, we provide nnmerical evidence showing the non- 
integrability of systems ffT^ - ffT3|) (corresponding to SR-geodesic flows 
with D = 6, 7, 8) by evolving the eqnation of motion of the rednced 
system ffTTj) . In Section |8] we have already claimed that system flT^ 
resembles the dynamics of a driven pendulnm that is chaotic. However, 
this resemblance can be a mere coincidence. In order to put forward a 
thorough investigation of whether in the above systems chaos appears 
or not, we need a more standardized method. 

One of the most classical methods for hnding chaos is given by in¬ 
vestigating the dynamics of the retnrn map on the snrfaces of section 
(Poincare map). We compnte this nnmerically by evolving the eqna- 
tions of motion with the Cash-Karp-Runge-Kutta scheme. The accn- 
racy of the nnmerical resnlts is checked by redncing the integration 
step size by an order of magnitnde and testing whether this rednction 
changes the trajectory of the orbit. 

Surfaces of section were employed in for the D = 8 case as 
wel@. There the snrface z = 0, z' > 0 has been chosen as the Poincare 


^Equations (80)-(82) of [Sa] with q = 0 correspond to our (fTTll with Q = Q 2 ix, y). 
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Figure 2. The surface of section on z = 0, 2 ' > 0 (black 
dots) for system ffT2]) with parameters a = 10, 6 = —0.1 
and the orbit starting from x = y = The red curve 
Qi{x,y) = 0 shows the limit of this section. The left 
panel shows 10® sections, while the right panel zooms 
the left plot around the starting point. 


section. However, if we employ the same surface of section for the 
D = 6 parabolic case we encounter a problem. Namely, the surface of 
section z = 0, z' > 0 does not meet all the trajectories in the phase 
space, because of the parabolic form of Qi(x, y). For example, in Fig. [2] 
the red curve Qi{x, 1 /) = 0 sets a limit for the section we can plot, and 
creates an obstacle to study the whole phase space. In other words, 
the surface z = 0, z' > 0 for system f[T2]) is not a good choice for the 
Poincare section. Moreover, the oscillations across the x-axis indicate 
that the system is non-compact. Namely, as an orbit evolves it tends 
to reach larger and larger values of |i/| and |a;|. 

Because of the above mentioned oscillations across the x-axis shown 
in Fig. [21 we assumed that a good surface of section for system flT2l) 
would be the surface x = 0, x' > 0. This assumption has proven 
to be correct and we show the results on Fig. [31 In both panels of 
Fig. [31 we can see a region of concentric closed curves (black curves), 
which represent regular orbits. The center of these regular orbits lies 
around the point {z, y) = (0, 20) in the left panel, and around the point 
{z, y) = {0, —20) in the right panel. The concentric curves indicate that 
the central point corresponds to a stable periodic orbit. 

In both cases around these concentric orbits he an irregular orbit 
(red dots), which tends to cover all the available phase space in the 
complement to the regular orbits. The irregular orbit apparently stems 
from a point around (z, y) = (0,30) in the left panel, and around 
{z^y) = (0,-30) in the right panel. Both these points match the 
appearance of unstable periodic orbits. 
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Figure 3. Details of the surface of section on x = 0, 
x' > 0 for system ffT^ with parameters a = 10 and 
b = —0.1 (left panel), 6 = 0.1 (right panel). The black 
closed curves represent regular orbits, while the red dots 
correspond to one irregular orbit. 


These irregular orbits in Fig. |3] indicate that the D = 6 parabolic 
system is non-integrable. Note that the plots of Fig. [3] do not show 
the whole phase space, because the system is non-compact. Instead we 
focus our plots on the region around the regular orbits and near the 
unstable point, where the irregular features are more prominent. 

System ffT3|) can be separated in two categories: the elliptic ones 
{ab > 0), and the hyperbolic ones {ab < 0). In the elliptic case the 
surface of section in Fig. 0] tells straightforwardly that the system is 
non-integrable. Namely, in Fig. 0] we can discern the characteristic 
features of a non-integrable system like chaotic regions and islands of 
stability belonging to Birkhoff chains. In the particular case of the 
system corresponding to D = 8, the indicated non-integrability is in 
agreement with the non-integrability conjecture of |Sa] . 

We can assert non-integrability also for the hyperbolic case on the 
ground of analytic dependence on the parameters a, b of our system 
(assuming the integrals should share the same property). However, 
we can conhrm this numerically as well, and we do it in Fig. 0) the 
hyperbolic orbit is shown on two different surfaces of section, and both 
of these surfaces indicate that the orbit is irregular, and therefore, the 
systems (in. corresponding to D = 7, 8, are non-integrable. 
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Figure 4. Dynamics on the snrface of section z = 0, 
z' > 0 for system with parameters a = 2, 6 = 1, 
c = 0 (elliptic case). 




Figure 5. The snrface of section on z = 0, z' > 0 (left 
and middle panel) for system flT^ with parameters a = 

— 10^, b = 10^. The orbit starts from x = y = 0, and 
evolves for 10^ sections. The whole orbit is shown in the 
left panel, while a detail focnsing on the starting point 
is shown in the middle panel. The right panel shows the 
same orbit on the snrface y = 0, y' > 0. 

Appendix A. SR-structures on 3D Lie groups 

Every left-invariant SR strnctnre on a 3-dimensional Lie gronp G is 
determined by a 2-dimensional snbspace (not snbalgebra) of the Lie 
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algebra 0 and a metric on it. The classification of such is due to [VGj . 
and this reference also contains the integration of the equations of 
geodesics in terms of a semi-direct product. 

Liouville integrability of left-invariant SR structures on 3D Lie groups 
G was proven in the preprint arXiv:math/0105128 of |Krj . It was later 
re-visited in |MS] . We provide a short proof here for completeness. 

Theorem 4. Non-holonomic geodesic flows of left-invariant SR-metrics 
on 3-dim Lie groups are Liouville integrable with polynomial integrals. 

Proof. The left-invariant Hamiltonian 2H = cjf -h CJ 2 commutes with 
all right-invariant forms di. Every 3-dimensional Lie algebra g has a 
Casimir function C G G°°(g*) (because G has odd dimension), so the 
involutive set of integrals is: I 2 = G, I 3 = 0i, where the number i is 
chosen such that R = H,l 2,13 are functionally independent. 

Moreover this G is linear for the Heisenberg algebra and quadratic for 
simple Lie algebras s[(2), so(3), but it can be non-algebraic (depending 
on parameters) in the remaining semi-direct cases g = xR^. In these 
cases, 1) = R^ is an Abelian subalgebra. The right-invariant forms R, I 3 
associated to a basis in f) are integrals in involution. The Hamiltonian 
H is algebraically (and functionally) independent of those, because 
otherwise it would be bi-invariant. This completes the proof. □ 

Appendix B. Prolongation of Killing symmetries 

Let g = © • • • © g_i be a (finite-dimensional) graded nilpotent 

Lie algebrqlj, such that g_i generates g. The Tanaka prolongation is 
a graded Lie algebra g such that g_ = ©*<001 = 0 and it is the maxi¬ 
mal graded Lie algebra with this property (its construction is outlined 
below). In particular, go = Oero( 0 ) is the algebra of grading preserving 
derivations of the Lie algebra g. 

Given a subalgebra go C go, the Tanaka prolongation pr(g, go) = 
Q-u © • • • © 0-1 © 00 © 01 © • • • is naturally defined if we restrict 
the non-positive part to g © go. Constructively gi consists of the 
homomorphisms ip : g_j —?■ gi_*, i > 0, satisfying the Leibniz rule 
<p([a;, y]) = [(p(x), i/] + [x, (p{y)], then we similarly define 02 etc. If some 
0 ,. = 0, then also g* = 0 for i > r and the algebra g is finite-dimensional. 

An example of reduction of go is given by a left-invariant SR-structure. 

Theorem 5. Let g be a Riemannian metric on g-i and go = Oero(0) H 
50(0-1,5-). Then pr_^(g, 0 o) = 0, i.e. g* = 0 Vi > 0. 

is customary in Tanaka theory to use negative gradation in the basic part, 
so we switch here from the notations used in the main body of the paper. 
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This theorem is due to Morimoto [Mo] . His proof is based on a result 
due to Yatsui. In the case of our interest we give a simpler argument. 

Proof in the case dimg_i = 2. Clearly the only possibility for 
non-zero go is M = 5o(2) = (cq) that acts on g_i = as a 

complex structure: [eo,e'_^] = e"i, [eo,e"i] = —ehi- 

For 0 7 ^ 9 ? G gi there is a basis of g_i such that = eo, </3(e" i) = 

0 . Then for e _2 = [e'_^,e"i] we have </5(e_2) = —ehi- Let e _2 = e _2 
and dehne recursively e_s = [e"]^,ei_s], s > 2. We have 9?(e_s) = ei_s 
so by induction e_s 7 ^ 0 Vs > 2, implying that dimg = 00 . This is a 
contradiction. □ 
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